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Abstract 

We study the relation between the quark distribution function q(x) and the frag- 
mentation function D q (z) based on a general form D q (x) = C(z)z a q(z) for valence 
and sea quarks. By adopting two known parametrizations of quark distributions for 
the proton, we find three simple options for the fragmentation functions that can pro- 
vide a good description of the available experimental data on proton production in 
e + e~ inelastic annihilation. These three options support the revised Gribov-Lipatov 
relation D q (z) = zq(z) at z — > 1, as an approximate relation for the connection be- 
tween distribution and fragmentation functions. The three options differ in the sea 
contributions and lead to distinct predictions for antiproton production in the reaction 
p + p — > p + X, thus they are distinguishable in future experiments at RHIC-BNL. 
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The quark distribution function q(x) and the quark fragmentation function D q (z) 
are two basic quantities on the structure of hadrons. It would be very useful if there 
exist simple connections between them, so that one can predict the poorly known 
D q (z) from the rather well known q(x). Conversely, one could predict the quark 
distribution functions of a hadron that cannot be used as a target from the quark 
fragmentation to the same hadron. In fact, some efforts have been made to connect 
fragmentation functions with the corresponding quark distributions [l], 0, || . There is 
a so called Gribov-Lipatov "reciprocity" relation |l| which connects the distribution 
and fragmentation functions in a form 

D q (z) oc q(x). (1) 

Such relation has been used as a useful "Ansatz" to model the quark fragmentation 
functions based on predictions of quark distributions functions |J . Recently, based 
on theoretical arguments with some assumptions @, a revised form of the Gribov- 
Lipatov "reciprocity" relation, i.e., 

D q (z) = zq(z), (2) 

has been suggested as an approximate relation at z — > 1. Thus it is necessary to check 
the validity of Eq. ([!]) and/or Eq. (g) by means of careful phenomeno logical studies. 

The nucleon is a satisfactory laboratory to check the relation between fragmenta- 
tion and distribution functions, since we have data both on the quark distributions of 
the nucleon from deep inelastic scattering (DIS) || || and on the fragmentation func- 
tions of quark to proton from e + e~ inelastic annihilation (IA). Various parametriza- 
tions of parton distributions have been obtained from the DIS experimental data 
on the nucleon, with rather high precision J7L |], [H|. We can start with any set 
of parton distributions to parametrize the fragmentation functions. In this paper, 
we will adopt the CTEQ parametrization (CTEQ6 set 1) [0] of the parton distribu- 
tions and for comparison, we will also use another recent parametrization obtained 
in the statistical physical picture, by Bourrely, Soffer and Buccella (BSB) 0. These 
parametrizations provide reliable information on the valence and sea quarks and their 
respective roles can be studied separately. On the other hand, some experimental data 
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on proton production in e + e IA are available ]TT| and can be used to constrain the 



shape of the fragmentation functions of quark to proton. By parametrizing the frag- 
mentation functions of quark to proton based on a reliable set of parton distributions 
and confronting them with the available experimental data on proton production, we 
thus have an effective way to learn about the relation between fragmentation and 
distribution functions. For this analysis, we adopt a general form to relate fragmen- 
tation and distribution functions and we make a distinction between valence and sea 
quarks, as follows 



The above forms are always correct, since C v (z) and C s (z) are in principle arbitrary. 
Strictly speaking, there is no way to discriminate between "valence" and "sea" frag- 
mentation, because both "valence" and "sea" quarks q (not q) can fragment at a same 
z. Since it is not possible to distinguish whether the fragmenting quark is "valence" 
or "sea", we should consider Eq. (Q) as a phenomenological parametrization for the 
fragmentation functions of quarks and antiquarks, as follows 



It will be shown in this Letter that, in order to fit the experimental data of proton 
production in IA, with fragmentation functions based on the above parametrizations, 
we can have very simple forms of C v (z) and C s (z) together with simple values for 
a. Three options are found to be very good: (1) C v = 1 and C s = for a = 0, (2) 
C v = C s = 1 for a = 0.5, and (3) C v — 1 and C s = 3 for a — 1. Since the sea quark 
contributions are negligible at large x, all of the three options support the revised 
Gribov-Lipatov "reciprocity" relation D q (z) = zq(z) at z — > 1, as an approximate 
relation for the connection between distribution and fragmentation functions. Option 
1 corresponds to the case of a suppressed sea, whereas Option 3 corresponds to the case 
of an enhanced sea, so they provide two extreme situations. Thus these three options 
lead to different fragmentation pictures related to the respective roles of valence and 
sea quarks. We will show that they provide different predictions for proton/antiproton 



D v (z) = C v (z)z a q v (z), 
D s {z) = C s (z)z a q s (z). 



(3) 



D q (z) = D v (z)+D s (z), 
D^z) = D s (z). 



(4) 
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productions in the reactions p + p — > p(p) + X, thus they are distinguishable in future 
experiments at RHIC-BNL. 

The revised Gribov-Lipatov relation Eq. (§J) is only known to be valid near 2 — > 1 
and on a certain energy scale Qq, in the leading approximation. Since distribution 
and fragmentation functions have different evolution behaviors, it is appropriate to 
apply the relation at a starting energy scale Qq, and then to evolve the fragmentation 
functions to the experimental scale. Our goal is to get a good description for the cross 
section of proton production in e + e~ IA [0]]. In our analysis, we adopt the initial 



energy scale Qq of the quark distributions, i.e., it is Qo = 1.3 GeV for the CTEQ 
parametrization J7j and Qq = 2 GeV for the BSB parametrization || . In addition, in 
the evolution the gluon fragmentation functions are also needed, and we assume that 
they have the same relation with the corresponding gluon distribution functions as 
the sea quarks. For the evolution we adopt the evolution package of Ref. ||, modified 
for the fragmentation functions. 

In the quark-parton model, the differential cross section for the semi-inclusive 
proton production process e + e _ — > p + X can be expressed to leading order as 

1 , EC q [D q (x E ,Q 2 ) + D q (x E ,Q 2 )] 

1 da „ ^ (5) 



ctot dx E J2 C q 

1 



where x E = 2E p /yfs, s is the total center-of-mass (cm.) energy squared, E p is the 
energy of the produced proton in the e + e~ cm. frame, and atat is the total cross 
section for the process. In Eq.(|J), C q reads 

C q = e 2 q - 2xiv e v q e q + X2(a 2 e + v 2 e )(a 2 q + v 2 ), (6) 



with 



1 s(s-M 2 y 



X2 



' 16 sin 2 6 W cos 2 6 W (s - M\) 2 + M 2 Z Y\ 

1 s 2 

256 sin 4 9 W cos 4 9 W (s - M|) 2 + M 2 Z Y\ ' 



(7) 



a e = -1, (9) 
-l+4sin 2 fl ty , (10) 
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a q = 2T 3q , (11) 

and 

v q = 2T 3q -4e q sm 2 9 w , (12) 

where T 3q = 1/2 for u, while T 3q = — 1/2 for d, s quarks. Moreover N c = 3 is the color 
number, e q is the charge of the quark in units of the proton charge, 9 is the angle 
between the outgoing quark and the incoming electron, 9w is the Weinberg angle, 
and Mz and V z are the mass and width of Z°. 

In Fig. |l| we present our results of the proton production in e + e~ IA with the 
three simple options for the relation between distribution and fragmentation functions 
Eq. ©. Option 1: C v = 1 and C s = for a = 0; Option 2: C v = C s = 1 for 
a = 0.5; and Option 3: C v = 1 and C s = 3 for a = 1. We find that all of these 
three options fit rather well the experimental data in a large z range, though there 
is some discrepancy at small z. From Fig. 1(a) and (b), we find that there are only 
small differences in the predictions when we use the two different choices of quark 
distributions parametrizations. Thus our analysis of the relation between distribution 
and fragmentation functions is not sensitive to the available parametrizations of quark 
distributions which are well constrained by a vast number of experimental data. All 
of the three options support the validity of the revised Gribov-Lipatov relation Eq. (fj) 
at z — > 1. 

As indicated above, the three options differ in the roles played by the valence and 
sea quark contributions in the fragmentation functions. Option 2 corresponds to a 
situation with equal contributions from valence and sea quarks, whereas Option 1 and 
Option 3 correspond to the situations with sea suppressed and enhanced respectively. 
Though such differences do not show up significantly in e + e~ — >• p + X, they will show 
up in the proton/antiproton productions in p + p — > p(p) + X, as will be shown in the 
following. Thus we can test these three different options by examining predictions 
in these processes, where the roles played by valence and sea quarks are different in 
comparison with e + e~ IA. For the sake of simplicity, we only adopt the fragmentation 
functions based on the CTEQ parametrization of the parton distributions. In fact, the 
new facility RHIC at BNL is expected to be running at some high cm. energy such 
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Figure 1: Predictions for the cross section of proton production in e + e~ annihilation. 
Dotted, solid, dashed curves correspond to the three options of 1, 2 and 3 respectively, 
with (a) CTQE and (b) BSB parametrizations of the parton distributions. The 
experimental data are taken from [[11], [12 
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as \fs = 200 GeV and yfs = 500 GeV. We thus use these fragmentation functions to 
predict the cross section of the proton/antiproton productions via 

p + p^H + X, (13) 

where the produced hadron H, is a proton or an antiproton, respectively. The invari- 
ant cross section of the above process can be expressed as 

e b-P- = T, Cdxj'dxMx^Q^hix^Q^D^z^Q^ — ^iab^cd). (14) 

" PH aMd JWa J *b dt 

In Fig. |2], we present the rapidity distribution of the cross section for proton produc- 
tion at the cm. energy = 200 GeV (left figure (a)) and = 500 GeV (right 
figure (b) ), for various values of the transverse momentum p? of the produced hadron. 
In Fig. |3|, the px distribution of the cross sections is shown at the above mentioned 
two values of the cm. energy. We see that the predictions do not show large differ- 
ences for the three options. In Fig. ^-^|, the rapidity and the transverse momentum 
distributions of the cross sections are shown for antiproton production and now we 
observe significant differences for the predictions from the three options. This can be 
easily understood since antiproton production is sensitive to the sea quark fragmen- 
tation. We thus conclude that antiproton production in the reaction p + p — > p + X 
can discriminate clearly between the three options of quark fragmentation functions. 

In summary, we have studied the relation between fragmentation and distribution 
functions based on known parametrizations of parton distributions and the available 
experimental data on proton production in e + e~ IA. With general relations between 
fragmentation and distribution functions, we find three simple options that can pro- 
vide a good description of the proton production data in e + e _ IA. All of the three 
options support the revised Gribov-Lipatov relation D q (z) = zq(z) at z — > 1 as an 
approximate relation, and they differ in the roles played by valence and sea quarks in 
the fragmentation. Such difference can show up especially from antiproton produc- 
tion in the reaction p + p — > p + X, thus they can be tested in future experiments at 
RHIC-BNL. 
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Figure 2: Predictions of the rapidity distribution of the cross section of proton pro- 
duction in pp collision. Dotted, solid, and dashed curves correspond to the predictions 
from the three options 1, 2 and 3 respectively, for cm. energy y/s = 200 GeV (a) 
and y/s = 500 GeV (b). 
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Figure 3: Predictions of the transverse momentum distribution of the cross section 
of proton production in pp collision. Dotted, solid, and dashed curves correspond 
to the predictions from the three options 1, 2 and 3 respectively, for cm. energy 
js = 200 GeV (a) and y/s = 500 GeV (b). 
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